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Abstract 

Processes leading to a final state with at least two high transverse momentum leptons 
(electrons or muons) are studied using the full e^p data sample collected by the HI exper- 
iment at HERA. The data correspond to an integrated luminosity of 463 pb^^. Di-lepton 
and tri-lepton event classes are investigated. Cross sections of the production of e+e~ and 
pairs are derived in a restricted phase space dominated by photon-photon collisions. 
In general, good agreement is found with Standard Model predictions. Events are observed 
with a total scalar sum of lepton transverse momenta above 100 GeV where the Stan- 
dard Model expectation is low. In this region, combining di-lepton and tri-lepton classes, 
five events are observed in e+p collisions, compared to a Standard Model expectation of 
0.96 ± 0.12, while no such event is observed in e~p data for 0.64 it 0.09 expected. 
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1 Introduction 



Within the Standard Model (SM) the production of multi-lepton events in electror0-proton 
collisions proceeds mainly via photon-photon interactions, 77 — i'+f^, where photons are 
radiated from the incident beam particles [1]. The clean experimental signature and the pre- 
cise SM prediction of this process provide a high sensitivity to searches for new phenomena 
producing multi-lepton final states. For example, the production of a doubly charged Higgs 
boson [2, 3] or processes involving generic bosons carrying two units of lepton number (bilep- 
tons) [4], could lead to multi-lepton events of large invariant mass. Measurements of both multi- 
electron [5] and muon pair [6] production at high transverse momentum (Pt) have already been 
performed by the HI Collaboration using a data sample corresponding to an integrated lumi- 
nosity of 115 pb~^. Events with high invariant mass M12 of the two highest Pt electrons were 
measured. Three events with two electrons and three events with three electrons were observed 
in the region M12 > 100 GeV, where the SM prediction is low [5]. 

In the present paper a measurement of multi-lepton production at high P^ using the complete 
e^p HERA collider data of the HI experiment is presented. The data are recorded at electron 
and proton beam energies of 27.6 GeV and 820 GeV or 920 GeV, corresponding to centre-of- 
mass energies ^/s of 301 GeV or 319 GeV, respectively. The total integrated luminosity of 
the data is 463 pb^^, which represents a factor of four increase with respect to the previous 
published results. The data comprise 178 pb^^ recorded in e^p collisions and 285 pb^^ in e+p 
collisions, of which 35 pb~^ were recorded at ^/s = 301 GeV. While previous measurements 
were performed mainly using e+p data, an e~p data set with more than a ten-fold increase in 
integrated luminosity is now analysed. The current analysis is extended with respect to those 
presented in [5, 6], such that all event final states with at least two high P^ leptons, electrons 
(e) or muons (yu), are now investigated. In addition, differential cross sections of the production 
of e~^e~ and /i"*"/!" pairs are measured in a restricted phase space dominated by photon-photon 
collisions. 



2 Standard Model Processes and their Simulation 

Multi-lepton events are simulated with the GRAPE [7] generator, which includes all elec- 
troweak matrix elements at tree level. The production mechanisms include 77, 7Z, ZZ in- 
teractions, internal photon conversion and the decay of virtual or real Z bosons. Different 
approaches for the cross section calculation are followed depending on the virtuality of the 
photon coupled to the proton and on the mass W of the hadronic final state. In the elastic region, 
— > e i'^i~p, the proton vertex is described in terms of dipole form factors. The quasi-elastic 
domain is defined hy W < 5 GeV or < 1 GeV^. In the region W < 2 GeV, a resonance 
parametrisation [8] is used for the proton vertex. In the remaining quasi-elastic phase space, 
a fit to photoproduction and deep inelastic scattering data is used [9]. In the inelastic regime, 
corresponding to electron-quark interactions with W > 5 GeV and > 1 GeV^, proton par- 
ton densities taken from the CTEQ5L [10] parametrisation are used. Initial and final state QED 

'in this paper the term "electron" is used generically to refer to both electrons and positrons, if not otherwise 
stated. 
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radiation is included. The hadronic final state is simulated via interfaces to PYTHIA [11] and 
SOPHIA [12] for the inelastic and quasi-elastic regimes, respectively. GRAPE predicts cross 
sections for ep — > e /i+/i^X and ep ^ e e+e"X processes, leading to eiifi and eee final states. 
Events with only two leptons (/^./x, e/x or ee) are observed if the scattered electron or one lepton 
of the pair is not detected. The ep — > e t'^t~X process with subsequent leptonic tau decays is 
also simulated with GRAPE and its contribution to the studied final states is found to be at most 
4%. 

Experimental backgrounds to ee and eee final states arise from events in which, in addition 
to a genuine electron, one or two fake electrons are reconstructed. Neutral current (NC) deep- 
inelastic scattering (DIS) events (ep eX) in which hadrons or radiated photons are wrongly 
identified as electrons constitute the dominant background contribution. QED Compton scatter- 
ing ep — > e^X may also contribute if the photon is misidentified as an electron. Background to 
the ejjL final state may arise from NC DIS events if hadrons are misidentified as muons. The NC 
DIS and Compton processes are simulated using the RAPGAP [13] and WABGEN [14] Monte 
Carlo (MC) generators, respectively. 

Generated events are passed through the full GEANT [15] based simulation of the HI appa- 
ratus, which takes into account the running conditions of the different data taking periods, and 
are reconstructed and analysed with the same program chain as is used for the data. 



3 Experimental Conditions 

A detailed description of the HI experiment can be found in [16]. Only the detector components 
relevant to the present analysis are briefly described here. The origin of the HI coordinate sys- 
tem is the nominal ep interaction point, with the direction of the proton beam defining the posi- 
tive z-axis (forward region). Transverse momentum is measured in the xy plane. The pseudora- 
pidity 7] is related to the polar angle 6'by?7 = — In tan(6'/2). The Liquid Argon (LAr) calorime- 
ter [17] covers the polar angle range 4° < 9 < 154° with full azimuthal acceptance. Electromag- 
netic shower energies are measured with a precision of a{E)/E = 12%/y/E/GeV ®1% and 
hadronic energies witha{E)/E = 50%/y/E/GeV®2%, as measured in test beams [18,19]. In 
the backward region, energy measurements are provided by a lead/scintillating-fiber (SpaCal) 
calorimeter [20] covering the range 155° < 6 < 178°. The central (20° < 9 < 160°) and for- 
ward (7° < 6* < 25°) inner tracking detectors are used to measure charged particle trajectories, 
to reconstruct the interaction vertex and to complement the measurement of hadronic energies. 
The LAr calorimeter and inner tracking detectors are enclosed in a super-conducting magnetic 
coil with a field strength of 1.16 T. From the curvature of charged particle trajectories in the 
magnetic field, the central tracking system provides transverse momentum measurements with 
a resolution of ap^/Px = 0.005PT/GeV © 0.015 [21]. The return yoke of the magnetic coil is 
the outermost part of the detector and is equipped with streamer tubes forming the central muon 
detector (4° < ^ < 171°). 

The luminosity is determined from the rate of the Bethe-Heitler process ep — > ep7, mea- 
sured using a photon detector located close to the beam pipe at z — — 103 m, in the backward 
direction. 



5 



Events having an electromagnetic deposit in the LAr with an energy greater than 10 GeV are 
detected by the LAr trigger with an efficiency close to 100% [22]. The muon trigger is based on 
single muon signatures from the central muon detector, combined with signals from the central 
tracking detector. The trigger efficiency is about 90% for di-muon events with muon transverse 
momenta larger than 10 and 5 GeV. A combination of the LAr and muon triggers is used if an 
electron and a muon are both present in the final state, resulting in an efficiency close to 100% 
for electron and muon transverse momenta above 5 GeV. 

4 Data Analysis 

4.1 Particle Identification 

Electron candidates are identified in the polar angle range 5° < ^ < 175° as a compact and 
isolated electromagnetic shower in either the LAr or SpaCal calorimeter. The electron energy 
and angular direction are measured by the calorimeters. The electron energy threshold is 5 GeV 
in the polar angle range 20° < 9 < 175° and is raised to 10 GeV in the forward region (5° < 
9 < 20°). The calorimetric energy measured within a distance in the pseudorapidity-azimuth 
{rj — (j)) plane R = \J ISjf + A^^ < 0.75 around the electron is required to be below 2.5% of 
the electron energy. In the region of angular overlap between the LAr and the central tracking 
detectors (20° < 9 < 150°), hereafter referred to as the central region, the calorimetric electron 
identification is complemented by tracking information. In this region it is required that a 
well measured track geometrically matches the centre-of-gravity of the electromagnetic cluster 
within a distance of closest approach of 12 cm. The electron is required to be isolated from 
any other well measured track by a distance R > 0.5 to the electron direction. Furthermore, 
the distance from the first measured track point in the central drift chambers to the beam axis 
is required to be below 30 cm in order to reject photons that convert late in the central tracker 
material. In addition, in the central region the transverse momentum of the associated electron 
track P|.** is required to match the calorimetric measurement Pf. such that 1 / — 1 /P|. < 
0.02 GeV~^ in order to reject hadronic showers. Due to a lower track reconstruction efficiency 
and higher showering probability, no track conditions are required for electron candidates in 
the forward (5° < 6* < 20°) and backward (150° < 9 < 175°) regions. In these regions, 
no distinction between electrons and photons is made. The resulting electron identification 
efficiency is 80% in the central region and larger than 95% in the forward and backward regions. 

Muon candidates are identified in the polar angle range 20° < 9 < 160° with a minimum 
transverse momentum of 2 GeV. Muon identification is based on the measurement of a track in 
the inner tracking systems associated with a track segment or an energy deposit in the central 
muon detector [6,23]. The muon momentum is measured from the inner track curvature in the 
solenoidal magnetic field. A muon candidate should have no more than 3.5 GeV deposited in a 
cylinder, centred on the muon track direction, of radius 25 cm and 50 cm in the electromagnetic 
and hadronic sections of the LAr calorimeter, respectively. Misidentified hadrons are strongly 
suppressed by requiring that the muon candidate be separated from any jet and from any track 
by P > 1. The efficiency to identify muons is ~ 90%. 

Calorimeter energy deposits and tracks not previously identified as electron or muon can- 
didates are used to form combined cluster-track objects, from which the hadronic energy is 
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reconstructed [24,25]. Jets are reconstructed from these combined cluster-track objects using 
an inclusive algorithm [26, 27] with a minimum transverse momentum of 2.5 GeV. 

4.2 Event Selection 

In order to remove background events induced by cosmic rays and other non-ep sources, the 
event vertex is required to be reconstructed within 35 cm in z of the nominal interaction point. 
In addition, to remove non-ep background topological filters and timing vetoes are applied. 

Multi-lepton events are selected by requiring at least two central (20° < 9 < 150°) elec- 
tron or muon candidates, of which one must have > 10 GeV and the other Pj,>5 GeV. 
Additional lepton candidates are identified in the detector according to the criteria defined in 
section HTTl All lepton candidates are required to be isolated with respect to each other by a min- 
imum distance in pseudorapidity-azimuth ofR> 0.5. Lepton candidates are ordered according 
to decreasing transverse momentum, P^' > P^'^^. Final states with all possible combinations 
of lepton candidates are investigated. Selected events are classified into independent samples 
according to the flavour and the number of lepton candidates (e.g. ee, e/x, efj.fi). 

In order to measure the lepton pair production cross section in a well defined region of phase 
space, sub-samples of ee and fxfx events dominated by photon-photon collisions are selected, 
labelled (77) e and (77) respectively. In these subsamples the two leptons are required to be 
of opposite charge and a significant deficit compared to the initial state must be observed in 
the differencell E — of the energy and the longitudinal momentum of all visible particles, 
E — Pz < 45 GeV. These two conditions ensure that the incident electron is lost in the beam 
pipe after radiating a quasi-real photon of squared four-momentum Q"^ lower than 1 GeV^. 

4.3 Systematic Uncertainties 

The following experimental systematic uncertainties are considered: 

• The uncertainty on the electromagnetic energy scale varies depending on the polar angle 
from 0.7% in the backward and central region to 2% in the forward region. The polar 
angle measurement uncertainty of electromagnetic clusters is 3 mrad. The identification 
efficiency of electrons is known with an uncertainty of 3 to 5%, depending on the polar 
angle. 

• The scale uncertainty on the transverse momentum of high Pt muons is 2.5%. The un- 
certainty on the reconstruction of the muon polar angle is 3 mrad. The identification 
efficiency of muons is known with an uncertainty of 5%. 

• The hadronic energy scale is known within 2% at high transverse momentum and 5% for 
events with a total hadronic transverse momentum below 8 GeV. 

^For fully contained events or events where only longitudinal momentum along the proton direction {+z) is 
undetected, one expects E — — 2E^ ~ 55.2 GeV, where E^^ is the energy of the incident electron. If the 
scattered electron is undetected, the threshold E — < 45 GeV corresponds to a cut on the fractional energy loss 
y^[E- P,)/2E° < 0.82. 
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• The uncertainty on the trigger efficiency is estimated to be 3% if at least one electron 
candidate is detected, and 6% if only muons are present in the final state. 

• The luminosity measurement has an uncertainty of 3%. 

The effect of the above systematic uncertainties on the SM expectation is determined by 
varying the experimental quantities by ±1 standard deviation in the MC samples and propagat- 
ing these variations through the whole analysis chain. 

Additional model systematic uncertainties are attributed to the SM Monte Carlo generators 
described in section [21 The theoretical uncertainty on the lepton pair production cross section 
calculated with GRAPE is 3%. The uncertainty on the QED Compton and NC DIS back- 
ground is 20%, as deduced from dedicated studies [5]. The total error on the SM prediction 
is determined by adding the effects of all model and experimental systematic uncertainties in 
quadrature. 



5 Results 

5.1 Multi-Lepton Event Samples 

The observed event yields are summarised in table [T] Only classes for which at least one data 
event is selected are shown. In all other classes, no event is observed and the SM prediction is 
also negligible. The observed numbers of events are in good agreement with the SM expecta- 
tions. The eee and ee channels are dominated by electron pair production. The e/i/i, /i/i and e/i 
channels contain mainly events from muon pair production. The e/i channel is populated if the 
scattered electron and only one muon of the pair is selected. Four data events are classified as 
ee/i compared to a SM expectation of 1.43 ± 0.26, dominated by the production of muon pair 
events where one muon is lost and the second electron candidate is due to a radiated photon. 
One event with four electron candidates is observed compared to a prediction of 0.33 ± 0.07. 
According to MC simulations, this signature is due to tri-electron events with a radiated photon. 

The distributions of the invariant mass M12 of the two highest leptons for the eee and 
e/i/i samples are shown in figures [Ua) and (b), respectively. The distributions of the invariant 
mass M12 of the two leptons in the di-lepton event classes are presented in figures [TJc), (d) and 
(e). An overall agreement with the SM prediction is observed in all cases. 

High invariant mass events (M12 > 100 GeV) are observed in the data. The corresponding 
observed and predicted event yields are summarised for all channels in table [2l The three ee, 
three eee and one /i/i high mass events have already been discussed extensively in previous HI 
publications [5,6]. No additional events in these classes are observed in the new data. One e/i 
and two e/i/i high mass events are observed in the new data. The e/i/i event with the largest 
invariant mass M12 = 127 ± 10 GeV is presented in figure [21 In this event, M12 is formed by 
the electron and the highest P-r muon. In the other e/i/i event, M12 is formed by the two muons. 
No eeee or eefi event is observed with a di-lepton invariant mass above 100 GeV, in agreement 
with the corresponding SM expectations below 0.01. 
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The results for e+p and e~p data are shown separately in table [2l The di-lepton (ee, efi and 
/i/i) and tri-lepton (eee and e/i/i) events at high mass M12 > 100 GeV are all observed in e+p 
collisions whereas no such event is observed in the e~p data. 

The topology of the multi-lepton events can be further investigated using the scalar sum of 
the lepton transverse momenta ^ P^. This variable indicates the "hardness" of the event and 
also offers a good sensitivity for new physics searches [28]. Figure [3] presents the correlation 
between ^ Pt and the invariant mass M12, separately for di-lepton and tri-lepton events. It can 
be seen that tri-lepton events may have a large M12 but only intermediate J2 ^t- In such cases, 
the high mass is formed by one forward and one backward lepton. However, for di-lepton event 
classes high M12 values also imply a large ^ Pt. 

The quantity ^ Pt can be used to select the most energetic events and allows the combina- 
tion of different topologies of di-lepton and tri-lepton events with electrons and muons. Figure|4] 
presents the distributions of ^ Px of the observed multi-lepton events compared to the SM ex- 
pectation. A good overall agreement between the data and the SM prediction is observed. For 
J2 Pt > 100 GeV, five events are observed in total, compared to 1.60 ± 0.20 expected from 
the SM (see table [3]). These five events were all recorded in the e~^p data, for which the SM 
expectation is 0.96 ± 0.12. Furthermore, the events correspond to the three ee and the two e/i/i 
events observed with M12 > 100 GeV. 



5.2 Cross Section Measurements 



Cross sections of the production of electron and muon pairs from photon-photon collisions are 
measured using the selected (77)2 and (77);^ samples. The kinematic domain of the measure- 
ment is defined by 20° < 6^^'^ < 150°, P^' > 10 GeV, P^^ > 5 GeV, < 1 GeV^ and 
y < 0.82. The data samples collected at y/s = 301 GeV and 319 GeV are combined taking 
into account their respective luminosities. Assuming a linear dependence of the cross section 
on the proton beam energy, as predicted by the SM, the resulting cross section corresponds 
to an effective ^/s = 318 GeV. The total numbers of observed (77)e and (77)^ events are in 
agreement with the SM expectations, as summarised in tablefT] In the (77) e sample, the contam- 
ination from NC DIS and QED Compton background events is 2%. No significant background 
is present in the (77) sample. 

The cross section is evaluated in each bin i using the formula 



where A^^'^^*'^ is the number of observed events in bin i, N^^"^ the expected contribution from 
background processes in bin i, C the integrated luminosity of the data and Ai the signal accep- 
tance in bin i. The signal acceptance is calculated using GRAPE MC events, as the ratio of the 
number of events reconstructed in bin i divided by the number of events generated in the same 
bin. It accounts for detection efficiencies and migrations between bins. The mean signal accep- 
tance is 45% for ep — i> e e'^'e^X events and 60% for ep e fi'^n^X events. The systematic 
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error of the measured cross section is determined by repeating the analysis after applying the 
appropriate variations to the MC for each systematic source, as described in section |431 

The measured ep ^ e e^e^X cross section integrated over the phase space defined above 
is (T = 0.67 ± 0.06 ± 0.05 pb, where the first error is statistical and the second systematic. 
The measured cross section for muon pair production, ep ^ e fi^ n^X, in the same phase 
space is a = 0.63 ± 0.05 ± 0.09 pb. The results are in agreement with the SM expectation 
of 0.63 ± 0.02 pb calculated using the GRAPE generator. Combining the di-electron and di- 
muon samples, an average lepton pair production cross section of cr = 0.65 ± 0.04 ± 0.06 pb is 
measured. 

The differential cross sections of lepton pair production as a function of the transverse mo- 
mentum of the leading lepton P!j^ , the invariant mass of the lepton pair M« and the hadronic 
transverse momentum are listed for each sample in table |4] and shown in figure [5] for the 
combined electron and muon samples. The measurements are in good agreement with the SM 
cross sections. 



6 Conclusion 

The production of multi-lepton (electron or muon) events at high transverse momenta is studied 
in e~^p and e~p scattering. The full e^p data sample collected by the HI experiment at HERA 
with an integrated luminosity of 463 pb^^ is analysed. The yields of di-lepton and tri-lepton 
events are in good agreement with the SM predictions. In each sample distributions of the 
invariant mass M12 of the two highest Pt leptons and of the scalar sum of the lepton transverse 
momenta ^ Pt are studied and found to be in good overall agreement with the SM expectation. 

Events are observed in the di-lepton and tri-lepton channels with high invariant masses 
M12 above 100 GeV. All such events are observed in e'^p collisions. Five of them have a 
J2Pt > 100 GeV, whereas the corresponding SM expectation for e~^p collisions is 0.96 ± 0.12. 

Differential cross sections for electron and muon pair production are measured in a restricted 
phase space dominated by photon-photon interactions. The measured cross sections are in 
agreement with the SM expectations. 
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Multi-Leptons at HERA (463 pb"^) 



Selection 


Data 


SM 


Pair Production (GRAPE) 


NC DIS + Compton 


ee 


368 


390 ± 46 


332 ± 26 


58 ±30 




201 


211 ±32 


211 ± 32 


< 0.005 




132 


128 ±9 


118 ±8 


10.0 ±2.5 


eee 


73 


70 ±7 


69.8 ± 7.0 


0.2 ±0.1 


e/iii 


97 


102 ± 14 


102 ± 14 


< 0.005 


ee.fj, 


4 


1.43 ±0.26 


1.18 ±0.20 


0.25 ±0.14 


eeee 


1 


0.33 ±0.07 


0.33 ±0.07 


< 0.005 


(77) e 


146 


138 ± 12 


135 ± 11 


3.0 ± 1.0 


(77)/. 


163 


162 ± 24 


162 ± 24 


< 0.005 



Table 1: Observed and predicted event yields for the different event classes. The errors on the 
predictions include model uncertainties and experimental systematic errors added in quadra- 
ture. The limits on the background estimations correspond to the selection of no event in the 
simulated sample and are quoted at 95% confidence level. 

Multi-Leptons at HERA (463 pb"^) 



M12 >100GeV 



Selection 


Data 


SM 


Pair Production (GRAPE) 


NC DIS -1- Compton 


All data (463 pb"^) 


ee 


3 


1.34 ±0.20 


0.83 ±0.11 


0.51 ±0.13 




1 


0.17 ±0.07 


0.17 ±0.07 


< 0.005 


e/i 


1 


0.59 ±0.06 


0.59 ±0.06 


< 0.005 


eee 


3 


0.66 ±0.09 


0.66 ±0.09 


< 0.005 


e/i/i 


2 


0.16 ±0.05 


0.16 ±0.05 


< 0.005 


e~^p collisions (285 pb ^) 


ee 


3 


0.76 ±0.11 


0.49 ±0.07 


0.27 ±0.07 




1 


0.10 ±0.04 


0.10 ±0.04 


< 0.005 


e/i 


1 


0.35 ± 0.04 


0.35 ± 0.04 


< 0.005 


eee 


3 


0.39 ±0.05 


0.39 ±0.05 


< 0.005 


e/i/i 


2 


0.09 ±0.03 


0.09 ±0.03 


< 0.005 


e p collisions (178 pb ^) 


ee 





0.58 ±0.09 


0.34 ± 0.04 


0.24 ±0.07 







0.07 ±0.03 


0.07 ±0.03 


< 0.005 


e/i 





0.24 ±0.03 


0.24 ±0.03 


< 0.005 


eee 





0.27 ±0.04 


0.27 ±0.04 


< 0.005 


e/i/i 





0.07 ±0.03 


0.07 ±0.03 


< 0.005 



Table 2: Observed and predicted multi-lepton event yields for masses M12 > 100 GeV for the 
different event classes in all analysed samples. The errors on the predictions include model 
uncertainties and experimental systematic errors added in quadrature. The limits on the back- 
ground estimations correspond to the selection of no event in the simulated sample and are 
quoted at 95% confidence level. 
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Multi-Leptons at HERA (463 pb"^) 


^Pt >100 GeV 


Data sample 


Data 


SM Pair Production (GRAPE) 


NC DIS + Compton 


e+p (285 pb-i) 


5 


0.96 ±0.12 0.78 ±0.09 


0.18 ±0.05 


e-p (178 pb"^) 





0.64 ±0.09 0.51 ±0.07 


0.13 ±0.04 


All (463 pb-i) 


5 


1.60 ±0.20 1.29 ±0.15 


0.31 ± 0.09 



Table 3: Observed and predicted multi-lepton event yields for J^Pt > 100 GeV. Di-lepton and 
tri-lepton events are combined. The errors on the predictions include model uncertainties and 
experimental systematic errors added in quadrature. 



Multi-Leptons at HERA (463 pb'O 



Variable 


Measured 


Measured 


Measured 


Pair Production 


range 


(e+e-) 


(/i+/i-) 


(average) 


(GRAPE) 


[GeV] 


[fb/GeV] 


[fb/GeV] 


[fb/GeV] 


[fb/GeV] 


rj, 




da/dP^^ 






[10, 15] 


108 ±10 ±8 


95 ±9 ±14 


101 ±7 ±9 


92 ±3 


[15,20] 


16 ±4 ±2 


23 ± 4 ±4 


20 ±3 ±2 


22 ±1 


[20, 25] 


7.4 ±2.5 ±0.6 


3.7 ±1.7 ±0.5 


5.5 ±1.5 ±0.5 


6.6 ±0.2 


[25, 50] 


0.60 ±0.30 ±0.06 


0.70 ±0.30 ±0.10 


0.70 ±0.20 ±0.07 


0.86 ±0.03 


Mu 




da/dMu 






[15,25] 


34 ±4 ±3 


29 ±3 ±4 


32 ±3 ±3 


28 ±1 


[25,40] 


16 ±2 ±2 


17 ±2 ±3 


17 ±2 ±2 


18.0 ±0.5 


[40, 60] 


2.7 ±0.8 ±0.2 


1.7 ±0.6 ±0.3 


2.2 ±0.5 ±0.2 


2.9 ±0.1 


[60, 100] 


0.20 ±0.15 ±0.02 


0.40 ±0.20 ±0.07 


0.30 ±0.13 ±0.04 


0.24 ±0.01 






da/dP^ 






[0,5] 


95 ±9 ±8 


94 ± 8 ±14 


94 ± 6 ±9 


93 ±3 


[5, 12] 


25 ±6 ±2 


13 ±4 ±2 


18 ±3 ±2 


15.4 ±0.5 


[12,25] 


2.8 ±1.3 ±0.3 


4.2 ±1.2 ±0.6 


3.8 ±0.9 ±0.4 


3.8 ±0.1 


[25, 50] 


0.20 ±0.30 ±0.03 


0.50 ±0.30 ±0.09 


0.40 ±0.20 ±0.06 


0.20 ±0.01 



Table 4: Differential cross sections da/dPj.^, da/dMu and da/dP^ averaged for each quoted 
interval for the process ep — > el'^C'X in a restricted phase space dominated by the photon- 
photon process (see text for details). Cross sections are measured for e^e~ or pairs. The 
cross section obtained from the combination of e+e" and iJi^id,' channels is also presented. The 
first error is statistical and the second is systematic. Theoretical predictions, calculated with 
GRAPE, are shown in the last column. 
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Multi-Leptons at HERA (463 pb"^) 




[GeV] 



Multi-Leptons at HERA (463 pb"^) 




Mee [GeV] 



Multi-Leptons at HERA (463 pb~^) 




Me, [GeV] 



Multi-Leptons at HERA (463 pb'^) 




[GeV] 



Multi-Leptons at HERA (463 pb"^) 




Mnn [GeV] 



Figure 1 : The distribution of the invariant mass of the two highest Pt leptons for events classi- 
fied as eee (a), efifi (b) and ee (c), nfi (d) and e/x (e). The points correspond to the observed data 
events and the open histogram to the SM expectation. The total error on the SM expectation 
is given by the shaded band. The component of the SM expectation arising from lepton pair 
production is given by the hatched histogram. 



15 



H1 



= 2.8 GeV 




X - 61 GeV 



Figure 2: Display (side view) of the e/x/i event observed with the highest M12 invariant mass. 
Indicated are the reconstructed tracks and the energy depositions in the calorimeters. The beam 
positrons enter the detector from the left and the protons from the right. 
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Figure 3: Correlation of the invariant mass M12 with the scalar sum of the transverse momenta 
for di-lepton (a) and tri-lepton (b) events. The bold dots represent the data while the 
small points represent the pair production (GRAPE) prediction for a luminosity ~ 1000 times 
higher than that of data. 
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Multi-Leptons at HERA (e*p and e"p, 463 pb"^) 




E Pj [GeV] 



Multi-Leptons at HERA (e"^p, 285 pb"") 




L Pj [GeV] 



Multi-Leptons at HERA (e"p, 178 pb"') 




E Pj [GeV] 

Figure 4: The distribution of the scalar sum of the transverse momenta Pt for combined 
di-lepton and tri-lepton event samples for all data (a) as well as for e^p (b) and e~p (c). The 
points correspond to the observed data events and the open histogram to the SM expectation. 
The total error on the SM expectation is given by the shaded band. The component of the SM 
expectation arising from lepton pair production is given by the hatched histogram. 
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Multi-Leptons at HERA (463 pb'^) 



> 

0) 

O 10-^ 
n 

Q. 



Q. 
T3 

"6 



10- 





H1 




(a) 


i-m-i 

r ^ 




' — I — 



10 15 20 25 30 35 40 45 50 

[GeV] 



^10-^. 

0) 

O 



O 10-3 ^ 



10-^ 



I — •-^ 



1-4 " ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ 



■■■■■■I- 



■■■■■■I- 



H1 

(b) 



l1i 



20 30 40 50 60 70 80 90 100 



> 
o 

O 
n 

X I- 

Q. 

■a 

D 



10-^ r 



10-=^ z 



10' 





H1 




(c) 


■ l-i-l 




' J ' 








: 1 





5 10 15 20 25 30 35 40 45 50 

[GeV] 



Pfi> 10 GeV, p/2>5GeV 

20° < 0^1'^^ < 150° 
y < 0.82, Q2 < 1 GeV^ 

i-|h hi Data 

SM (GRAPE) 



Figure 5: The measured cross section for lepton pair production in a restricted phase space dom- 
inated by the photon-photon process as a function of the leading lepton transverse momentum 
P^^ (a), the invariant mass of the lepton pair (b) and the hadronic transverse momentum 
(c). The differential cross section is averaged over the intervals shown. The inner error bars 
represent the statistical errors, the outer error bars the statistical and systematic errors added in 
quadrature. The bands represent the SM prediction with its one standard deviation uncertainty. 
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